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Paramagnetic centers in a solid-state environment usually give rise to inhomogenously broadened 
electron paramagnetic resonance (EPR) lines, making conventionally detected free induction decay 
(FID) signals disappear within the spectrometer dead time. Here, experimental results of an elec- 
trically detected FID of phosphorus donors in silicon epilayers with natural isotope composition 
are presented, showing Ramsey fringes within the first 150 ns. An analytical model is developed to 
account for the data obtained as well as for the results of analogous two-pulse echo experiments. 
The results of a numerical calculation are further presented to assess the capability of the method 
to study spin-spin interactions. 



I. INTRODUCTION 

Solid-state based quantum computer (QC) architec- 
tures have recently attracted considerable interest due 
to the prospect of possible scalability. Several imple- 
mentations have been proposed, such as superconducting 
tunnel junctions [Ij, nuclear spins in crystal lattices [2j, 
nuclear spins of donors in Si |3j, and electron spins lo- 
calized in quantum dots (QDs) [iJ [5] or at donors [6j. 
However, a realistic solid-state environment introduces 
an intrinsic inhomogeneity of the qubit properties due 
to impurities, defects, interfaces etc., resulting in pos- 
sible computational errors. Nevertheless, fault-tolerant 
error correction schemes promise to compensate such er- 
rors up to a certain level [THS], making a quantification 
of the inhomogeneity a crucial step in the assessment of 
the suitability of specific physical realizations of qubits 
for QC. 

In the case of QDs and donors, the Zeeman levels of 
the electron spin in an external magnetic field are used 
as the qubit states. The Zeeman energy of each electron 
spin will vary within a spin ensemble, e.g. due to the ran- 
dom orientation of nuclear spins in the lattice that gen- 
erate a spatially fluctuating local magnetic field at the 
position of each electron spin. The resulting distribution 
^(cjs) of the Larmor frequencies cjs of the electron spins 
can be quantified in the frequency domain by the inho- 
mogenous linewidth obtained from a traditional contin- 
uous wave (cw) electron paramagnetic resonance (EPR) 
experiment |1Q|. Alternatively, a time domain analysis 
of ^(cjs) can be carried out by measuring the decay of 
the transverse magnetization in the free induction decay 
(FID) of a conventional pulsed EPR experiment. The 
problem encountered in both methods is their respective 
detection limit which is far above the number of spins 



* corresponding author, email: hoehne@wsi. turn. del 



present in nanoscale structures proposed for QCs. More- 
over, FID has not seen a widespread application for stud- 
ies of paramagnetic species in solids since the broad Lar- 
mor frequency distribution usually encountered in solid- 
state spin systems makes the FID signal disappear on a 
time scale that is usually much shorter than the spec- 
trometer dead time to after a microwave pulse of typi- 
cally 100 - 200 ns [UHIS]. 

Optically and electrically detected magnetic resonance 
experiments (ODMR and EDMR, respectively) have 
shown much promise since they exceed conventional 
EPR concerning detection sensitivity by at least six or- 
ders of magnitude [HI [15]. In the case of phospho- 
rus donors (^^P) in silicon, the development of pulsed 
EDMR (pEDMR) pHTS] has paved the way for study- 
ing dynamical properties of spin systems by applying 
different microwave pulse sequences, such as using the 
Hahn-echo sequence to study coherence times T2 [19], 
or the inversion recovery sequence to study spin-lattice 
relaxation times Ti [20j. The underlying principle of 
pEDMR on ^^P donors is the use of microwave pulses 
for the generation of coherent states of ^^P spins, which 
can be detected by transient photoconductivity measure- 
ments [17, 18j, making use of spin-dependent recombi- 
nation of spin pairs formed by ^^P and paramagnetic 
Si/Si02 interface states [211122] . 

In current pulsed EPR, the dead-time problem of FID 
has been circumvented by several techniques such as 
high-frequency EPR [13j, FID-detected hole burning [llj 
or detection of the EPR spectra by the Electron Spin 
Echo (ESE), although the latter is more sensitive to 
distortions through nuclear modulations in contrast to 
FID [IT, 13] . In this paper, we study the possibility of 
electrically detecting (ED) FID using a tomography tech- 
nique developed earlier [19] to investigate the Larmor fre- 
quency distribution ^(cjs) of a ^^P ensemble in natural 
silicon ("^*Si). It is shown that information within the 
usual dead time of conventional EPR-detected FID can 
be obtained. An analytical equation is deduced to de- 
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scribe the experimental data, which in turn agrees weh 
with the results of continuous wave (cw) EDMR experi- 
ments. Furthermore, a numerical study is performed to 
assess the EDFID technique in terms of its capability to 
quantitatively investigate the coupling of spin pairs. 



II. EXPERIMENTAL DETAILS 

The sample used in this work is fabricated by chemi- 
cal vapor deposition and consists of a 22 nm thick "^*Si 
layer with [P] = 9 x 10^^ cm~^ covered by a native oxide. 
It is grown on a 2.5 pm thick, nominally undoped "^*Si 
buffer on a silicon-on-insulator substrate. Evaporated in- 
terdigit Cr/Au contacts with a period of 20 pm covering 
an active area of 2 x 2.25 mm^ are biased with 100 mV, 
resulting in a current of ^ 80 pA under illumination with 
white light. The major paramagnetic states are identified 
in cwEDMR experiments performed in a Bruker X-band 
dielectric microwave resonator for pulsed EPR. The mi- 
crowave frequency of = a;mw/(27r) = 9.7400 GHz 
is provided by an HP83640 microwave synthesizer. The 
measurements are performed at 5.5 K in a helium gas fiow 
cryostat. The samples are oriented in an external mag- 
netic field ^0 with the [001] axis of the Si wafer parallel 
to ^o- The first derivative of the relative current change 
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Figure 1: First derivative spectrum of the relative change 
AI/I of the photocurrent in a cwEDMR experiment on 
^^*Si:P with a magnetic field modulation amplitude of 0.2 mT. 
The red curve represents a fit using Gaussian lineshapes and 
equal amplitudes of both "^^P resonance lines. The constituent 
lines of the fit are shown below for better visibility. 

AI/I [cf. Fig.[l] is measured with magnetic field modu- 
lation and lock-in detection as a function of Bq provided 
by a Bruker BE 25 electromagnet. We used the two pro- 
nounced ^^P hyperfine resonance lines (denoted hf(^^P)) 
at Bo = 346.17 mT and Bq = 350.27 mT to calibrate 
the magnetic field such that their center field corre- 
sponds to ^ = 1.9985 [23j. In addition, a resonance line 
at ^ = 2.0069 ± 0.0004 arising from the Pbo center at the 
Si/Si02 interface can be observed in accordance with pre- 



vious studies for Bq \ \ [001] [53]. A further resonance line is 
observed at ^ = 2. 0036 ±0.0004 which, due to its ^-factor, 
is attributed to the Pbi center [25j. We are aware of the 
fact that there is no concensus in the literature whether 
the Pbi defect is electrically active [26 , 27J and therefore, 
measurements of the angular dependence of the ^-factor 
or the hyperfine interactions with ^^Si nuclei would be 
needed for an unambiguous identification. The small 
resonance line at = 348.12 mT could be attributed 
to conduction band electrons (CE) with a ^-factor of 
g ^ 1.9990 ± 0.0004 [28J. A fit assuming Gaussian line 
shapes and equal amplitudes of both ^^P resonance lines 
is indicated by the red line in Fig.[T] The decomposition 
of the fitted spectrum is shown in the lower part of Fig.[T] 
An analysis of the isolated high-field hf(^^P) line yields a 
peak-to-peak line width of A5f,^™^^ = 0.30 ± 0.03 mT 
after correcting for the infiuence of magnetic field modu- 
lation [29j. This line broadening is predominantly inho- 
mogeneous, caused by randomly oriented nuclear spins 
of the ^^Si isotope in natural Si, which give rise to an 
unresolved superhyperfine multiplet [30j. Homogeneous 
broadening only plays a minor role since Ti and T2 times 
of the ^^P donor spins determined by ED inversion re- 
covery [20J and ED Hahn echo [19j experiments on this 
sample yield Ti ^ 5.3 ps and T2 ^ 3.3 ps, respectively, 
when the measured decays are fitted by a single expo- 
nential dependence. 

The pulsed EDMR experiments are performed at a mi- 
crowave frequency of i^mw = ^mw/ (27r) = 9.7331 GHz un- 
der the same orientation of the sample as the cwEDMR 
experiments. For all pulsed measurements, the magnetic 
field ^0 is corrected by the magnetic field offset deter- 
mined by cwEDMR. The microwave pulses are shaped 
using a SPINCORE PulseBlasterESR-Pro 400 MHz pulse 
generator and a system of microwave mixers, and are 
then amplified by an Applied Systems Engineering 117X 
traveling wave tube with a maximum peak power of 1 kW. 
The actual pulse shapes coupled into the resonator are 
checked in refiection experiments. The quality factor of 
the dielectric resonator is adjusted to make a compro- 
mise between sufficient excitation bandwidth and a mi- 
crowave magnetic field Bi high enough for coherent spin 
manipulation. The adjustment of the microwave power is 
achieved by a tuneable attenuator; the 7r/2-pulse time of 
^7r/2 = 15 ns corresponding to Bi = 0.6 mT used through- 
out this paper is determined in Rabi-oscillation experi- 
ments as previously developed [ITl [18] . The current tran- 
sients are acquired using a current amplifier followed by 
a voltage amplifier and a fast data aquisition card. To 
improve the signal-to- noise ratio, we applied a two-step 
phase cycling sequence where the phase of the last 7r/2 
pulse was switched by 180° at a frequency of ~ 5 Hz and 
the signals for each phase are then subtracted from each 
other. Switching the phases at a frequency of several Hz 
has the same purpose as a lock-in detection scheme and 
improves the signal-to-noise ratio by reducing the effects 
of 1 //-noise in the measurement setup. 
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III. RESULTS AND DISCUSSION 



A. Electrically detected FID 
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Figure 2: Electrically detected free induction decay or Ram- 
sey experiment on the high- field hf('^^P) resonance, (a) Con- 
tour plot of AQ as a function of the external magnetic field 
Bo and the free evolution time r. White dashed curves mark 
the positions of local extrema of AQ described by Eq. 
The black dashed line indicates the position where the cross- 
sectional chart shown in panel (b) is taken, (b) Cross-section 
of AQ along the evolution time axis at Bo — 349.25 mT. Red 
dashed curves illustrate the Gaussian-shape damping of the 
oscillation amplitude. 

The EDFID tomography is performed by a 7r/2-r-7r/2 
pulse sequence with varying evolution time r, consisting 
of the conventional free induction pulse sequence 7r/2-r 
followed by a 7r/2-projection pulse as usually applied in 
multi-pulse EDMR experiments [19, 20j. Hence, it coin- 
cides with the pulse sequence of the Ramsey experiment 
[31j. The photocurrent transients measured for each Bq 
and r are integrated over a fixed recording time inter- 
val satisfying the criteria described in Ref. [32j , resulting 
in a charge difference AQ. Figure|2] shows experimen- 
tal results of an EDFID tomography experiment on the 
isolated high- field hf(^^P) line. 

We will now show that the pattern in Fig.[2](a), which 
is characteristic of an EDFID can be understood by a 
simple model in which the contribution of the state of 
each spin pair at the end of the second 7r/2 pulse is 
proportional to its projection onto the singlet state l^*) 



[T71 [T5] . Hence, the measured charge Q cx —S^^{r) = 
— Tr(|5')(5'|p) reveals the average singlet content of the 
spin pair ensemble described by the density operator p. 
This is in contrast to conventional ESR, where for an 
FID the magnetization after a 7r/2 pulse is detected. For 
microwave frequencies close to the Larmor frequency of 
the high-field hf(^-'^P), the singlet content S{r) of each 
spin pair reflects the dynamics of only this spin species 
[TT] while in a first approximation the spin state of Pbo 
is unaltered and just serves as a projection partner. This 
is justified since the separation of the Larmor frequencies 
of the ^^P and Pbo spins for the high- field hf(^^P) res- 
onance is approximately one order of magnitude larger 
than the on-resonance Rabi frequency uji = gji^Bi/h. 
The minor effects of the off-resonance excitation of the 
other resonance lines can be seen as small oscillations on 
the Ramsey pattern in Fig.[2](a) at magnetic fields lower 
than 350 mT. We also neglect spin-spin interaction and 
incoherent processes during the pulse sequence. The for- 
mer will be addressed in Sec. lIII Cl and the latter is a valid 
assumption since the time constant for the fastest inco- 
herent process is T2 ~ 3.3 ps as measured in ED Hahn 
echo decay experiments on this sample. With these as- 
sumptions, an expression for the theoretically expected 
signal 
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with uji = cji/v^ and Auj = cjq — ^mw can be derived fol- 
lowing Ref. [33, 34J as shown in AppendixjAj In Eq. Q, 
quantifies the width of th e Larmor frequency distri- 
bution as defined in Eq. (A5). The locations of the local 
extrema of AQ are given by Eq. (|A8|) as 
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representing hyperbolas in the Bo-r-pldiue. These hy- 
perbolas fit the experimentally observed pattern well, 
as evident from the white dashed curves in Fig. [2] (a). 
The exponential term in Eq. ([T]) describes an envelope 
in the time domain which shows an exp [— (r/TFID)^]- 
type decay behavior [35j as depicted by red dashed line 
in Fig. [2] (b) with the time constant 
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This implies that for short pulses, i.e. in the high 
microwave power limit, TpiD is inversely proportional 
to the width of the Larmor frequency distribution 
and thus AB^^^^^^ of the M{^^F) resonance line in 
cwEDMR. The actual decay characteristics deviate from 
the exp [— (t/Tfid)^] behavior since the lineshape is 
a convolution of a Gaussian and Lorentzian lineshape 
rather than a pure Gaussian. 

In Fig.[3](a), cross-sections of AQ along the evolution 
time axis taken at different values of Bq are plotted as 
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Figure 3: Oscillations in the evolution time domain, (a) 
Cross-sections of AQ along the evolution time axis taken at 
different values of Bq. The damping of the oscillations indi- 
cates that the spin ensemble dephases. Red curves are fits 
based on the model given by Eq. Q . (b) Frequency z^fid of 
the damped oscillations as a function of the external magnetic 
field Bo shows a linear behavior as expected from Eq. 



a function of the evolution time r, revealing strongly 
damped oscillations. These characteristics are consistent 
with those described by Eq. ([T]) . A clear dependence of 
the frequencies of the damped oscillations z^fid on the ex- 
ternal magnetic field Bq can be observed. The red lines 
in Fig.jsj^a) show fits of the oscillations by the function 



AQ ~ —Ae V^fidJ cos [27ri/FiD(T + tq)] 



(4) 



which is based on the model given in Eq. ([T]). In all the 
fits shown in Fig.|3](a), a global phase correction of tq ^ 
20 ns has to be taken into account which is comparable 
to the overall pulse length of 30 ns. It can be attributed 
to the fact that dephasing during the pulse times can 
not be neglected due to the finite pulse width compared 
to the evolution time r. The values of vfib obtained 
from the fits are plotted as a function of the external 
magnetic field Bq and displayed in Fig.[3](b). A clear 
linear dependence of z^fid on the magnetic field Bq can 
be observed as expected from Eq. ([T]), 
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with k = gjiiB/[h{l-h2{a^/uJi)^)]. This is consistent with 
the linear dependence of the oscillation frequency on the 
detuning in a Ramsey experiment [361 EZ]- From a lin- 
ear fit of the data the value of 5res = 350.02 ±0.01 mT 
is obtained, which corresponds to the center position of 
the high-field hf(^^P) resonance line. Using Eq. ^ and 
27r/a;i = 60 ns, the average value Tfid = 64.4 ± 5 ns 
obtained from the fits of the damped oscillations can 
be related to an expected cwEDMR peak-to-peak line 
width of AB™ = 2ha^/{g/j.B) = 0.26 ± 0.02 mT. This 
is in agreement with the result obtained from cwEDMR 
A5f,^™^^ = 0.30 ±0.03 mT, demonstrating the consis- 
tency of the experiments. 



B. Electrically detected Hahn echo.jpg 

Electrically detected echo sequences have been previ- 
ously used to study T2 [19] and Ti times [20j of phospho- 
rus donors near Si/Si02 interface defects. In this section, 
we will present detailed experimental results of the elec- 
trically detected Hahn echo measurement with a focus on 
the fine structure of the echo response. 

The echo is measured using the previously developed 
tomography technique [19j extending the pulse sequence 
7r/2-ri-7r-r2 of the conventional two-pulse spin echo con- 
taining two free evolution times ri and r2 by a final 7r/2 
pulse as applied in the EDFID technique. The measure- 
ments are conducted under the same experimental con- 
ditions as the EDFID experiments. 
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Figure 4: Electrically detected two-pulse Hahn echo on the 
high- field hf(^^P) resonance with n = 300 ns. (a) Contour 
plot of AQ as a function of the external magnetic field Bo 
and the free evolution time T2. Black dashed lines mark 
the hyperbola pattern according to Eq. ([6|. Red and blue 
dashed lines indicate positions where cross-sectional diagrams 
shown in panel (b) are taken. (b) Cross-section of AQ 
along the evolution time axis at resonance condition (red 
line, Bo — 350.0 mT) and for the off resonant case (blue line, 
Bo — 350.5 mT). Please refer to the text for details. 

Figure|4] shows experimental results of an ED spin echo 
tomography experiment on the isolated high- field hf(^^P) 
line with t\ = 300 ns held fix. The values of AQ plotted 
as a function of B^ and T2 [Fig.jlja)] are obtained in the 
same way as described in the EDFID section. Cross- 
sections along the evolution time axis are displayed in 
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Fig.gb). The red curve taken at the resonance field 
Bq = 350.0 mT shows a Gaussian-shaped peak centered 
around r2 = ri = 300 ns. The cross section of AQ at the 
off-resonance field Bq = 350.5 mT, which is represented 
by the blue curve, shows oscillations as a function of the 
evolution time r2 with a maximum at r2 = ri = 300 ns, 
which decay for |r2 — ti| ^ 100 ns. 

The characteristic pattern indicated by the black 
dashed hyperbolas can be understood by the same quan- 
titative model described in the previous section on ED- 
FID. The singlet content S{r) proportional to the re- 
combination probability P-^^i of a single spin after a 7r/2- 
Ti-7r-T2-7r/2 pulse sequence can be calculated using the 
matrix formalism described in Ref. [33l [34]. For a Lar- 
mor frequency distribution modelled by a Gaussian with 
standard deviation <Jc^ centered about cjq, we can derive 



analogously to Eq. (A7) that 
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with uji = uJi/2 
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Figure 6: (a) Experimental results of the EDFID mea- 
surement showing all spectral features, (b) Simulation for 
J = MHz. In both simulation and experimental results, 
the characteristic Ramsey pattern can be clearly seen at the 
position of the high- field hf(^^P) resonance, whereas the pat- 
terns at the low- field hf(^"'^P) and the Pbo resonances are more 
complicated due to mutual interference. 
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Figure 5: Oscillations in the evolution time domain, (a) 
Cross-sections of AQ along the evolution time axis taken at 
different values of Bq. Red curves are fits based on the model 
given by Eq. ([6|. (b) The frequency z^Echo of the damped 
oscillations as a function of the external magnetic field Bq 
shows a linear behavior as expected from Eq. ([7|. 



cwEDMR linewidth of A5ff,^° 
value is consistent with the values A5pp^^^^ 
0.03 mT and AB™ 



From the fits of the data, an average value of Wecho = 
106.9 ± 4 ns is obtained, corresponding to an expected 

0.28 ± 0.01 mT. This 
0.30 ± 

0.26 ± 0.02 mT obtained from 
previous experiments within the accuracy limits. There- 
fore, the ED Hahn echo response on the high- field hf(^^P) 
resonance shows the same fine structure as the EDFID 
experiment and can be explained by the same model as 
expected from the fact that the echo pulse sequence con- 
sists of two FIDs back to back [39j. However, the small 
oscillations seen in EDFID are not observed in the ED 
Hahn echo response, since they are fully defocussed due 
to the additional central 7r-pulse. 



C. Spin-spin coupling 



Similar to the analysis of the EDFID experiment, var- 
ious cross-sections of AQ along the evolution time axis 
are taken and shown in Fig.jsja). Data fitting based on 
Eq. (|6| is performed and illustrated by the red curves. 
The oscillation frequency i^Echo obtained from the fits is 
plotted as a function of the external magnetic field Bq 
[Fig.jsjb)], where a linear dependence of i^Echo as a func- 
tion of ^0 can be observed as expected from 

1 o'^ 

Z^Echo 7^ I -2 (^0 ~ ^mw) = k' {Bq - ^res) , (7) 

with = gii^/[h{l + 4(crc^/co'i)^)] analogous to Eq. ([s]). 
The exponential term in Eq. (|6| describes a Gaussian en- 



So far the experimental results have been discussed 
in the context of off-resonance oscillations and dephas- 
ing due to inhomogeneous line broadening. In different 
previous studies [321 0^ , the possible impact of coupling 
between the partners of the spin pair on EDMR exper- 
iments has been discussed. In the following, results of 
a numerical study of the EDFID experiment discussed 
above are presented, focussing on the possibility of ED- 
FID to estimate the coupling strength. 

The system is modelled by an ensemble of spin S = 1/2 
pairs described by the density operator p. The Hamilto- 
nian of an individual pair is defined as 

n = no^nj^ni{t) (9) 
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= -^P/iB (^0 ± 5hf/2 + 5shf) 
+ 2^dbMB (^0 + ^Adb) erf" 



(10) 



representing the static uncoupled Hamiltonian in the 
presence of a constant magnetic field Bq = BqGz super- 
imposed with the hyperfine field of ^^P 5hf = 4.2 mT 
and the superhyperfine field 5shf at the position of 
the donor, where the latter can be considered fixed for 
timescales shorter than the precession period of ^^Si nu- 
cleus [41j. ^Adb is the local shift of the static magnetic 
field at the position of the Pbo center due to effects such 
as disorder and superhyperfine interactions. The cr^^y^z 
denote the Pauli spin operators. The circularly polarized 
microwave of angular frequency cjmw and magnitude Bi 
is represented in the rotating frame by 



(11) 



which is nonzero during the pulse. Spin-spin interaction 
is modelled by an exchange coupling Hamiltonian repre- 
sented by 
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with & = {cFxi^yi^z)^ ' The dipolar coupling is smaller 
than 1 MHz for interspin distances larger than 3 nm, 
which is neglected here for simplicity. The simulation of 
the spin pair ensemble dynamics is based on the Liouville 
equation dtp = i[p^i-L]~ /hin which, in contrast to Eq. (5) 
of Ref. [17J, all terms related to incoherent processes are 
dropped since the time constant of the fastest incoher- 
ent process is more than one order of magnitude larger 
than the duration of the pulse sequence, as already men- 
tioned in Section H. The initial steady state of the den- 
sity operator is assumed to be given by the pure triplet 
state pS = (|T+)(T+| + |T_)(T_|) /2 with |T+) = It^t^) 
and |T_) = H^^i-^) [II]. For triplet recombination rates 
tt much smaller than the singlet recombination rate rs, 
the observable Q{t) refiecting the state of the pair en- 
semble at the end of the second 7r/2 pulse assumes the 
form Q{t) oc —5{t) = — ((^p^db^p + (^p^db^p) [32j, where 
—^Pii = —{pii{^) ~ Pii) denotes the negative difference 
between the diagonal elements of the density matrix at 
the end of the second 7r/2-pulse and the initial steady 
state. The negative sign expresses the quenching of the 
photocurrent due to recombination. Inhomogeneous line 
broadening is taken into account by calculating —S{r) for 
a single spin pair and subsequent averaging over Gaus- 
sian distributions for both ^shf and ^Adb with exper- 
imentally obtained standard deviations from the pulsed 
EDFID spectrum shown in Fig.[6j Furthermore, the sim- 
ulation takes all four combinations of spin pair forma- 
tion into account, which arise from the two resonance 
positions of ^^P, the Pbo, and the Pbi- 

Figure[6]^a) shows the complete experimental results 
encompassing all resonances of the EDFID tomography 
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Figure 7: (a) and (b) Fourier transformation of the EDFID 
tomography experiment and simulation shown in Fig. [6] The 
linear dependence of the oscillation frequency on Bq — Bres 
described by Eq. ^ is clearly seen in the frequency domain. 
The two ^^P resonances and the broader Pbo resonance are 
marked by black and white dashed lines, respectively, (c) and 
(d) For better visibility, details of the FFT-spectrum near 
the high-field M{^^F) resonance are shown after subtracting a 
background of the form of Eq. Q from the experimental and 
simulated data. Again, the low-field hf(^^P) resonance and 
the Pbo resonance are marked by black and white dashed lines. 
An additional resonance, indicated by the solid white line, can 
be seen in the experimental data. Its spectral position is in 
accordance with the small central line observed in cwEDMR 
(see Fig.[T]). 



experiment discussed in Fig.|2ja). Figurejof^b) shows the 
simulation of —(5 for J = as a function of Bq and r af- 
ter subtraction of a constant background obtained from 
the value of —5 for large r, resulting in the quantity A 
which can be compared to AQ in the experiment. The 
characteristic patterns of simulation and experiment fit 
quite well. At the high-field ^^P resonance, small oscil- 
lations superimposed on the Ramsey oscillation pattern 
can be seen in the experimental data as well as in the sim- 
ulation. These small oscillations are due to the partial 
excitation of the low-field ^^P and the Pbo spins by the 
microwave pulses on the high-field ^^P resonance. Details 
of these patterns are shown in the Fourier transformed 
data shown in Fig.[7| (a) and (b). The linear dependence 
of the oscillation frequency on Bq — ^res described by 
Eq. ([5]) is clearly visible in the frequency domain. The 
two P resonances and the Pbo resonance are marked 
by black and white dashed lines, respectively. The Pbi 
resonance is not resolved due to the spectral overlap with 
the Pbo and its smaller amplitude. For better visibility, 
details of the FFT-spectrum near the high-field hf(^^P) 
resonance are shown in panels (c) and (d) after subtract- 
ing a background of the form of Eq. (|4| from the experi- 
mental and simulated data. Again, the low-field hf(^^P) 
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Figure 8: Simulation of the EDFID experiment at the high- 
field hf(^^P) resonance for different exchange coupling param- 
eter. An oscillation of the signal in the time-domain with the 
coupling frequency is expected. This oscillation is masked 
by the exponential decay of the signal due to dephasing and 
therefore can only be resolved for J > 5 MHz. 
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resonance and the Pbo resonance are marked by black and 
white dashed lines. An additional resonance, indicated 
by the solid white line, can be seen in the experimental 
data. Its spectral position is in accordance with the small 
central line observed in cwEDMR (see Fig.[T|. This line 
is not taken into account in the simulation. The intensity 
of the Fourier amplitude of the lines as a function of the 
magnetic field can b e described by an equation o f th e 
form sin^ (|Vl + j (1 + x^)), analogous to Eq. ( [ai| . 
In particular, the minima at frequencies of ~ 65 MHz 
and ~ 130 MHz correspond to rotations of the spins by 
integer multiples of 2i\. 

The simulation is further extended to nonzero coupling 
parameters with the focus on the clearly observable pat- 
tern structure on the high- field hf(^^P) resonance, which 
is shown in Fig. [8) Starting from J = 0, the exchange 
coupling parameter is increased in steps of 5 MHz result- 
ing in a change of the qualitative behavior of the char- 
acteristic pattern. Whereas for J = each hyperbola 
(cf. Fig.[2| either indicates positions of local maxima or 
minima, the values of A on each hyperbola oscillate as 
a function of r for J 7^ 0. This behavior can be clearly 
observed on the axis of symmetry at 5o = ^res, which is 
displayed in Fig.|9|a) for different values of J. 

For vanishing coupling, A relaxes exponentially to the 
equilibrium. For larger values of J, damped oscillations 
of A with frequency J are formed decaying to the equi- 
librium within the dephasing time TpiD- Compared with 
the existing experimental data [cf. Fig.[3](a)], the cou- 
pling is estimated <C 5 MHz. Please note, that in these 
simulations only the high-field ^^P and the Pbo are taken 



Figure 9: (a) Cross section along the r axis at the high- field 
resonance field of "^^P of the experimental data (open circles) 
and the simulated data (solid lines) for different exchange 
coupling parameter between ^^P and Pbo- (b) Simulated data 
like in (a) for an isotopically purified ^^Si sample. Here, more 
oscillations caused by the weak coupling can be seen due to 
the longer dephasing time. 



into account. Therefore, the small oscillations in Fig. [9] 
are due to off-resonance excitations of the Pbo spins. 

Clearer insight might be obtained in an experiment us- 
ing isotopically purified ^^Si samples as shown in the sim- 
ulation in Fig.[9]^b) where the line broadening (although 
determined by a homogeneous line width) is qualitatively 
modelled by a Gaussian distribution with FWHM line 
width of 0.023 mT corresponding to a concentration of 
29Si nuclei of ^ % [30j. Since TpiD ^l^u^ [cf. Eq. ^\ 
oscillations on the axis of symmetry could be observed 
within the dephasing time also for weak J. Qualitatively, 
the frequency of the oscillations increases with increasing 
J as evident from Fig.joj^b). In the general case, distri- 
butions of exchange interaction as well as dipolar inter- 
action have to be taken into account [17J, which would 
result in an averaging out of the oscillation pattern. How- 
ever, previous studies have revealed that only spin-pairs 
within a narrow range of intra-pair distances contribute 
to the observed signals [20j, which should make an esti- 
mate of J of contributing spin pairs still possible. 
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IV. SUMMARY 

To summarize, we have used pulsed EDMR to study 
the free induction decay of phosphorus donor spins in 
silicon. We can resolve oscillations up to 150 ns lim- 
ited by dephasing due to superhyperfine interactions with 
surrounding ^^Si nuclei. An analytical model is used to 
describe the FID of an inhomogeneously broadened line 
which is in good agreement with the experimental data. 
In addition, structures on two-pulse electron spin echoes 
have been measured which can be described by the same 
analytical model. The results of a numerical calculation 
are further presented and compared with the experimen- 
tal data to assess the capability of the method to study 
spin-spin interactions. From these results, we can give an 
upper bound for the coupling parameter of J ~ 5 MHz 
in the samples studied. 
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For distributions ^(cjs) with a maximum at the center 
frequency cJo, the dominant term of S^^{r) for Umw close 
to ujQ is given by 



sin^ (fVl + a;2) 1 + cos (Ar) 



(A3) 



with X = A/wi, neglecting a term of 0{\oJ mw — t^sP)- To 
obtain an analytical expression, Eq. ( A3 ) can be further 
simplified by the approximation 



sin^ (lyiT^) 



(A4) 



since both functions share the same leading orders in 
the Taylor expansion, tolerating a deviation of 6% in the 
integrated area within the interval defined by the zero- 
crossings of sin^(...)/(l + x^) in Eq. (A3). Modelling the 
Larmor frequency distribution by a Gaussian 



exp 



(A5) 



Appendix A: analytical expression describing the 
edfid pattern 



with standard deviation and center cjq, the average 
singlet content is given by 



In this section, Eq. ([2| used to describe the pattern in 
Fig.[2|a) is derived. Neglecting spin-spin interactions and 
incoherent processes, the singlet content S{t) is propor- 
tional to the fiipping probability P-^^^ of a single spin after 
a 7r/2-r-7r/2 pulse sequence, which has been investigated 
in studies related to nuclear magnetic resonance [33 , 34j 



P^^^(r) = 4sin^l9sin^ 



- cos t/ sm 



2 

Ar 
T 



sm 



/At 
cos(- 



atr 



COS 



atr 



(Al) 



with 



sm( 



a 



where cjs is the Larmor frequency of the ^^P donor elec- 
tron and cjmw the microwave frequency, = 7r/(2co'i) 
and r denote the length of the 7r/2 pulse and the free 
evolution time, respectively. For inhomogenously broad- 
ened lines, the observable 6'^^(r) is obtained by averaging 
S{t) (X i{r) over the Larmor frequency distribution 
[17J 





" 1 Aw2 - 




oc exp 




1 1 + exp 



X COS 



~,2 



(A6) 



with uji = cJi/V^ and Auj = ujq —uJmw Since Q cx —5^^, 
AQ is proportional to - [S^^'lr) - 6'^^(r oo)] as the 
constant background given by 5'^^(r ^ oo) is identical 
for the signals obtained for both phases (+x and -x) of the 
last 7r/2 pulse and thus subtracted by the data evaluation 
procedure described in Sec. [TT| This results in 



AQ cx — exp 



1 ^l^l 



2a^ 



r,2 



tAujt 



(yf. 



(A7) 



with local extrema approximately determined by values 
of Po and r for which the cosine term in Eq. ( A7 ) is equal 
to ±1, i.e. 



Po — 



n7rn(l + 2(aa;/cji)2) 1 



neZ. (A8) 



oc 

/ 



This term represents hyperbolas in the Bq-t plane shown 
(A2) in Fig.[2l 
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